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The  research  supported  by  this  grant  included  two  projects:  /.  Temperature  Control  In 
Polymer  Extrusion  Processes,  here  an  optimization  problem  was  formulated.  This  problem 
was  motivated  by  the  desire  to  obtain  uniform  extrudate  temperature  at  the  die  exit  in  a 
polymer  extrusion  process.  Control  was  effected  by  adjustments  to  the  heat  flux  along  the 
surface  of  the  pipe.  An  optimality  system  of  partial  differential  equations  was  derived  from 
which  optimal  controls  and  states  may  be  determined.  Then,  finite  element  discretizations 
of  the  optimality  system  were  defined  and  error  estimates  were  provided  along  with  an 
efficient  solution  algorithm  for  the  discreet.  Finally,  computational  results  were  given  for  a 
model  example  with  Oldroyd  type  fluid,  demonstrating  the  effectiveness  of  our  theory  and 
methods,  as  well  as  their  potential  applicability  to  industrial  problems.  II.  Analysis  and 
Finite  Element  Approximation  of  An  Optimal  Control  Problem  in  Electrochemistry  with 
Current  Density  Controls,  here  an  optimal  control  problem  for  impressed  cathodic  systems 
in  electrochemistry  was  studied.  The  control  in  this  problem  was  the  current  density  on  the 
anode.  A  matching  objective  functional  was  considered.  The  existence  of  an  optimal 
solution  was  proved.  The  use  of  Lagrange  multiplier  rules  was  justified  and  an  optimality 
system  of  equations  established.  Finally,  a  finite  element  algorithm  was  defined  and 
optimal  error  estimates  were  derived. 


Temperature  Control  In  Polymer  Extrusion  Processes 


In  polymer  extrusion  processes,  one  is  often  interested  in  maintaining  a  quasi-uniform 
temperature,  to  reduce  material  (extrusion  product)  inhomogeneity,  throughout  the 
extrudate  cooling  process.  In  this  project,  we  studied  the  somewhat  simplified  case:  one 
tries  to  obtain  a  uniform  temperature  distribution  at  the  exit  under  a  steady  state  situation. 
The  means  by  which  we  achieved  such  a  uniform  temperature  distribution  at  the  exit  were 
to  adjust  the  heat  flux  on  the  surface  of  the  pipe  near  the  exit. 

The  extrudate  in  equation  was  assumed  to  be  a  viscoelastic  fluid  of  the  Oldroyd  type  such 
as  polymer  melts  with  a  fast  relaxation  mode.  We  took  as  the  governing  equations  for  the 
Oldroyd  type  the  Navier-Stokes  equations,  the  incompressibility  constraint,  the  appropriate 
constitutive  equation  for  an  Oldroyd  model,  the  energy  equation,  and  simplified  boundary 
conditions.  The  control  function  was  determined  in  a  manner  that  would  allow  hot  spots  to 
be  avoided. 

We  investigated  two  means  of  obtaining  a  uniform  temperature  distribution.  The  first  is  to 
make  the  gradient  of  the  temperature  along  a  portion  of  the  boundary  small.  Another  means 
of  achieving  the  desired  result  is  to  try  to  force  the  temperature  field  itself  to  be  quasi¬ 
uniform.  Numerical  experiments  revealed  that  both  techniques  work  effectively  for  the 
desired  objective.  However,  we  focus  our  research  on  the  latter. 

We  proved  the  existence  and  uniqueness  of  optimal  solutions  and  derived  an  optimality 
system,  that  is,  a  set  of  equations  from  which  the  optimal  control  and  state  may  be 
determined.  Also,  finite  element  methods  and  numerical  examples  were  presented.  We  also 
developed  an  iterative  algorithm  to  compute  the  approximate  solution.  The  convergence  of 
our  algorithm  was  proved  and  a  comparison  with  a  direct  method  made. 


Analysis  and  Finite  Element  Approximation 
of  an  Optimal  Control  Problem  in 
Electrochemistry  with  Current  Density  Controls 

We  investigated  an  optimal  control  problem  for  impressed  cathodic  systems.  A  typical 
example  of  an  impressed  cathodic  system  is  a  metal  container  filled  with  an  electrolyte.  The 
painted  portion  of  the  container  surface  is  usually  treated  as  insulated.  The  unpainted  apart 
is  divided  into  cathode  and  anode  that  are  connected  to  the  negative  and  positive  poles  of  an 
electrical  source,  respectively.  By  adjusting  the  current  density  on  the  anode  we  could 
effectively  alter  the  potential  distribution  on  the  entire  bounding  surface  or  in  the  entire  flow 
domain.  The  potential  distribution,  of  course,  affects  on  the  chemical  reaction  process 
occurring  inside  the  flow  domain,  which  in  turn  affect  the  rate  of  corrosion  of  the  metal 
container.  Thus,  the  current  density  on  the  anode  can  be  used  as  a  practical  control  variable 
for  generation  a  desired  potential  field.  This  idea  can  be  conveniently  formulated  as  optimal 
control  problems  for  the  potential  equation  with  appropriate  boundary  conditions.  Optimal 
control  problems  of  this  sort  have  been  studied.  Here  existence  of  an  optimal  solution  is 
proved.  The  use  of  Lagrange  multiplier  rules  is  justified  and  an  optimality  system  of 
equations  is  established.  Finally,  a  finite  element  algorithm  is  defined  and  optimal  error 
estimates  are  derived. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


Abstract.  An  optimization  problem  is  formulated  motivated  by  the  desire 
to  obtain  uniform  extrudate  terrperature  at  the  die  exit  in  a  polymer 
extrusion  process.  Control  is  effected  by  adjustments  to  the  heat  flux 
along  the  surface  of  the  pipe.  An  optimality  system  of  partial  differential 
equations  is  derived  from  v\^ich  optimal  controls  and  states  may  be 
determined.  Then,  finite  element  discretizations  of  the  optimality  system 
are  defined  and  error  estimates  are  provided  along  with  an  efficient 
solution  algorithm  for  the  discrete  equations.  Finally,  corputational 
results  are  given  for  a  model  exaitple  with  Oldroyd  type  fluid,  demonstrating 
the  effectiveness  of  our  theory  and  method  as  well  as  their  potential 
applicability  to  industrial  probl^ns. 


BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


In  polymer  extrusion  processes,  one  is  often  interested  in  maintaining  a  quasi-uniform 
temperature,  so  as  to  reduce  material  (extrusion  product)  inhomegeneity,  tl^oughout  the 
extrudate  cooling  process.  In  this  research,  we  study  how  one  obtains  a  uniform 
temperature  distribution  at  the  exit  under  a  steady  state  situation.  The  means  we  use  to 
achieve  such  a  uniform  temperature  distribution  at  the  exit,  is  to  adjust  the  heat  flux  on  the 
surface  of  the  pipe  near  the  exit. 

The  extrude  in  question  is  assumed  to  be  a  viscoelastic  fluid  of  the  Oldroyd  type,  such  as 
polymer  melts  with  a  fast  relaxation  mode.  Two  means  of  obtaining  a  uniform  temperature 
distribution,  resulting  in  two  different  functionals,  have  been  investigated.  The  first 
involves  making  the  gradient  of  the  temperature  along  the  boundary  small.  Thus,  for 
example,  given  a  velocity  field ,  we  seek  a  temperature  field  and  a  control  field  such  that 
the  functional 


M(T,g)=|-  IgradTfdT  +  J^  IgfdT 

Jro  Jr^ 

is  minimized  subject,  of  course,  to  the  constraints  imposed  by  the  flow  equations.  Here  the 
minimization  results  in  a  quasi-uniform  temperature  distribution  along  the  boundary 
segment,  since  the  surface  derivatives  of  the  temperature  are  forced  to  be  small.  Another 
means  of  achieving  the  desired  result  is  to  try  to  erectly  force  the  temperature  field  itself  to 
be  quasi-uniform.  Thus,  now,  given  a  velocity  field,  we  would  seek  a  temperature  field 
and  a  control  field  such  that  the  second  functional 

N(T.g)=i[ 

Jro 

is  minimized  subject  to  the  constraints  imposed  by  the  flow  equations.  The  non-negative 
parameters  y  and  5  can  be  used  to  change  the  relative  importance  of  the  two  terms  appearing 
in  the  definition  of  N  as  well  as  to  act  as  penalty  parameters.  Our  numerical  experiments 
demonstrated  that  a  small  y  is  more  useful  in  achieving  quasi-uniform  boundary 
temperature  distributions,  although  it  also  reduces  the  accuracy  of  the  approximate  solution. 
Numerical  experiments  also  show  that  both  M  and  N  work  effectively  for  the  desired 
objective.  However,  our  work  has  focused  on  N.  Under  the  realistic  assumption  that  u  •  n 
=  o  on  the  boundary,  we  have  proved  the  existence  and  uniqueness  of  optimal  solutions 
and  derived  an  optimality  system,  i.e.,  a  set  of  equations  fi’om  which  the  optimal  control 
and  state  may  be  determined.  In  addition,  finite  element  methods  have  been  used  to 
compute  an  approximate  solution  of  the  optimality  system.  Optimal  error  estimates  have 
been  derived  and  numerical  experiments  have  been  performed.  Finally,  we  have  developed 
an  iterating  algorithm  to  compute  the  approximate  solution.  The  convergence  of  our 
algorithm  has  been  proved  and  a  comparison  with  the  direct  method  made. 

Control  problems  for  the  fully  coupled  problem  as  well  as  temperature  matching  for  the 
entire  extrudate  will  be  addressed  next 
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Abstract.  An  opliiiiization  problem  is  formulated  motivated  by  the  desire  to  obtain  unifonu  extrudate  lemiwraturu  at 
the  die  exit  in  a  polymer  extrusion  process.  Control  is  effected  by  adjustments  to  the  heat  flux  along  the  surface  of  (he  i>i|)e. 
An  optimality  system  of  partial  differential  equations  is  derived  from  wliich  optimal  controls  and  states  may  be  delerniined. 
Then,  finite  element  discretizations  of  the  optimality  system  are  defined  and  error  estimates  are  provided  along  with  an  efiicient 
solution  algorithm  fur  the  discrete  equations.  Finally,  computational  results  are  given  for  a  model  example  with  Olilroyd  type 
lluiil,  demonstrating  the  effectiveness  of  our  theory  and  method  as  well  as  their  potential  applicability  to  industrial  problems. 

1.  Iiitrotliictioii 

III  ])olymer  extrusion  processes,  one  is  often  interested  in  maintaining  a  quasi-uniform  temperature,  so 
as  to  reduce  material  (extrusion  product)  inhomogeneity,  throughout  the  extrudate  cooling  process.  In  this 
[taper,  we  will  study  the  somewhat  simplihed  case:  one  tries  to  obtain  a  uniform  temperature  distribution  at 
the  exit  under  steady  state  situation.  The  means  we  use  to  achieve  such  a  uniform  temperature  distribution 
at  the  exit  is  to  adjust  the  heat  flux  on  the  surface  of  the  pipe  near  the  exit. 

'I'lie  extrudate  in  question  will  be  assumed  to  be  viscoelastic  fluid  of  Oldroyd  type  such  as  polymer  melts 
with  a  fast  relaxation  mode  [10].  Let  u  denote  the  velocity  field,  p  the  pressure  field,  T  the  temperature  field 
and  T  the  purely  elastic  part  of  the  extra  stress.  Let  D  =  ^(gradii-l-  gradu^),  W  =  ^(gradu  —  grad  ) 
and  DaT  =  (u  grad  )t-1-tW  —  Wr  —  a(DT-l-TD)  where  —  1  <  a  <  1-  The  parameters  Re,  We  and  u>  are  the 
Reynolds  number,  Weissenberg  number  and  retardation  parameter,  respectively.  The  governing  equations 
for  the  Oldroyd  type  fluid,  in  dimensionless  form,  is  given  by  the  Navier-Stokes  equations 

ide(u  •  grad  )u -t- grad  p  =  (1  —  w)Au -1- divT-f  f ,  in  fl ,  (1.1) 

the  incompressiblity  constraint 

divu  =  0  in  n ,  (1.2) 

t  The  work  of  MDG  was  supported  by  the  Air  Force  Oflice  of  Scientific  Research  under  grant  number  AFOSR-90-0179  and 
by  the  Oilice  of  Naval  Research  under  grant  number  N000I4-91-J-1493i  that  of  LSI!  was  supported  by  the  Natural  Science  and 
Fngiiieei'ing  Research  Council  of  Canada  under  grant  number  OGP-0089763;  that  of  JCT  was  supported  by  US  Army  Research 
Oflice  under  grant  number  DAAL03-91-G0237. 


the  constitutive  equation  (Oldroyd  model) 


r  +  WeDaT=  2uD  infi, 


(1.3) 


and,  for  simplicity,  the  boundary  condition 


u  =  h  on  r , 


(1.4) 


and  also  the  energy  equation 

— K  AT  +  (ugrad  )T 

=  Q  +  2/i(grad  u  +  gradu^)  :  (grad  u  +  grad  u^)  in  Q , 


(1.5) 


with  boundary  conditions 


T  =  0  on  Ff) , 

^  =  Hn  onFArUFo, 
an 


(1.6) 

(1.7) 


—  =if  onFc-  (1-8) 

The  data  functions  f,  Q,  Hn,  and  h  are  assumed  known;  the  control  5  is  to  be  determined  so  that  hot  spots 
are  avoided.  The  constants  k  and  /i  depend  on  the  thermal  conductivity  coefficient,  density,  specific  heat 
at  constant  volume,  and  viscosity  coefficient  of  the  fluid.  See  [11]  for  details.  We  assumed  that  buoyancy 
effects  can  be  negelected,  and  thus  the  temperature  variable  does  not  appear  in  (1.1). 

Two  means  of  obtaining  a  uniform  temperature  distribution  come  to  mind.  The  first  is  to  make  the 
gradient  of  the  temperature  along  the  boundary  F,,  small.  Thus,  for  example,  given  a  velocity  field  u,  we 
would  seek  a  temperature  field  T  and  a  control  field  g  such  that  the  functional 


M{T,g)  =  ^J^  |grad,r|2dF+y^J<,pdF 


(1.9) 


is  minimized  subject,  of  course,  to  the  constraints  imposed  by  the  flow  equations  (1.5)-(1.8).  Here,  grad^ 
denotes  the  surface  gradient  operator,  e.g.,  in  the  tangential  derivative  operator  d/dr.  The  non-negative 
parameters  q  and  6  can  be  used  to  change  the  relative  importance  of  the  the  two  terms  appearing  in  the 
definition  of  .V(  as  well  as  to  act  as  penalty  parameters.  The  appearance  of  the  control  g  in  the  definition  of 
is  necessary  because  we  are  not  imposing  any  a  priori  limits  on  the  size  of  this  control.  The  minimization 
of  (1.9)  results  in  a  quasi-uniform  temperature  distribution  along  the  boundary  segment  F<,  because  the 
surface  derivatives  of  the  temperature  are  forced  to  be  smtill.  Another  means  of  achieving  the  desired  result 
is  to  try  to  directly  force  the  temperature  field  itself  to  be  quasi-uniform.  Thus,  now,  given  a  velocity  field 
u,  we  would  seek  a  temperature  field  T  and  a  control  field  g  such  that  the  functioned 


JiT,9)  =  ^f  \T-T,\^dr  +  !^  [  |<7pdF 

JTo  ^  JVc 


(1.10) 


is  minimized  subject  to  (1.1)-(1.8),  where  Td  is  some  desired  temperature  distribution,  e.g.,  something  close 
to  the  average  temperature  along  Fc  for  the  uncontrolled  system.  The  non-negative  parameters  7  and  6  can 
be  used  to  change  the  relative  importance  of  the  the  two  terms  appearing  in  the  definition  of  as  well  as  to 
act  as  penalty  parameters.  As  will  be  demonstrated  by  numerical  examples  in  §6,  a  small  7  is  more  useful 
in  achieving  quasi-uniform  boundary  temperature  distributions,  although  it  also  reduces  the  accuracy  of  the 


approximate  solution.  We  will  examine  the  latter  issue  in  §4.  Numerical  experiments  show  that  both  (1.9) 
and  (1.10)  work  effectively  for  the  desired  objective.  We  will  focus  our  discussion  on  (1.10)  thoroughout  this 
paper. 

Under  the  realistic  assumption  that  u  •  n  =  0  on  Tc  U  Tjv,  we  may  prove  the  existence  and  uniqueness 
of  optimal  solutions  and  derive  an  optimality  system,  i.e.,  a  set  of  equations  from  which  the  optimal  control 
and  state  may  be  determined.  Also,  finite  element  methods  are  used  to  compute  an  approximate  solution 
of  the  optimality  system.  Optimal  error  estimates  are  derived  and  numerical  examples  are  presented.  We 
have  also  developed  an  iterative  algorithm  to  compute  the  approximate  solution.  The  convergence  of  our 
algorithm  is  proved  and  a  comparison  with  the  direct  method  is  made. 

Control  problems  for  the  fully  coupled  problem  as  well  as  temperature  matching  for  the  entire  extrudate 
will  be  addressed  elsewhere. 

We  close  this  section  by  introducing  some  of  the  notation  used  in  subsequent  sections.  Throughout,  C 
will  denote  a  positive  constant  whose  meaning  and  value  changes  with  context.  Also,  H‘{V),  s  G  II ,  denotes 
the  standard  Sobolev  space  of  order  s  with  respect  to  the  set  D,  where  T)  is  either  the  flow  domain  Q,  or  its 
boundary  T,  or  part  of  that  boundary.  Of  course,  H°{V)  -  L'^{V).  Dual  spaces  will  be  denoted  by  (•)*.  Of 
particular  interest  will  be  the  space 

H]){Q.)  =  {S  e  H\Q.)  :  5  =  0  onTi,}. 

Norms  of  functions  belonging  to  H’{r)  and  /f*(rc)  are  denoted  by  ||  •  ||j  ,  ||  •  ||,,r  and  ||  •  ||i,rci 

respectively. 

The  inner  product  in  L^(Q)  is  denoted  by  (•,•),  that  in  i^(r)  by  (^^ri  that  in  L^{ro)  by  (-(  Oro) 
that  in  L^iTc)  by  (-.  Orc-  Since,  in  general,  we  will  use  T^-spaces  as  pivot  spaces,  these  notation  will  also 
be  employed  to  denote  pairings  between  Sobolev  spaces  and  their  duals. 

We  will  use  the  bilinear  form 

a{T,S)=  f  gradT  grad5dn  Vr,5ei/^(D) 

Jn 

and  the  trilinear  form 

c(u,r,5)  =  /(u-gradr)5dQ  VuGH^fi)  and 
Jn 

These  forms  are  continuous  in  the  sense  that  there  exist  constants  Cg  and  Cc  >  0  such  that 

|a(r,5)|<c„||T|ji||5||i  ^T,SeH\Q)  (1.11) 

and 

|c(u,r,5)|<c,||u||i||rili|1511i  VuGH‘(fi)  and  WT,SeH\n).  (1.12) 

Moreover,  we  have  the  coercivity  property 

a(r,r)  >  Caiiriif  vr g  (i.i3) 

for  some  constant  Cg  >  0. 

For  details  concerning  the  notation  employed  and  the  inequalities  (1.11)-(1.13),  one  may  consult,  e.g., 
[1]  and  [7]. 
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2.  The  Optimization  Problem,  Existence  of  Solutions,  And  Optimality  System 

We  begin  by  giving  a  precise  statement  of  the  optimization  problem  we  consider.  We  will  assume  the 
domain  fi  is  a  polygon  in  We  first  recall  that  (1.1)-(1.4)  uncouples  from  (1.5)-(1.8).  We  may  solve  for 
(u,pt)  from  (1.1)-(1.4)  once  and  for  all  and  then  plug  them  into  (1.5)-(1.8).  Thus  the  only  state  variable  is  T, 
i.e.,  the  temperature  field,  and  the  only  boundary  control  variable  is  g.  The  state  and  control  variables  are 
constrained  to  satisfy  the  system  (1.5)-(1.8),  which  we  recast  into  the  following  weak  form:  find  T  € 
such  that 

Ka{T,S)-¥c{vi,T,S)  =  {Q,S)  +  K{g,S)rc^>^{HN,S)rs  V5€/fi,(n),  (2.1) 

where  we  have  introduced  the  simplifying  notation 

Q  =  Q  +  2/i(grad  u  +  grad  u^)  :  (grad  u  +  grad  u^) . 

Note  that  since  we  seek  T  €  Hq{Q), 

T  =  0  on  Tp  .  (2.2) 

Throughout,  we  will  assume  that  the  given  velocity  field  u  is  smooth  and  satisfies 

divu  =  0  inf)  and  un>0  a.e.  on  Fp  U  Fat  .  (2-3) 

Under  these  assumptions,  we  have  the  useful  relation 

c(u,  5, 5)  =  ^  /  (u  •  n)S^  dF  >  0  ^  S  E  H\Q.)  ,  (2.4) 

2  7r 

which  can  be  derived  by  setting  T  =  S  in  the  following  integration  by  parts  formula: 

c(u,r,5)  =  J  iu-ii)TSdr-ciu,S,T). 

For  each  possible  control  function  g,  there  exists  a  unique  corresponding  state  function  T. 

Lemma  2.1-  For  every  g  €  I<^(Fc),  there  exists  a  unique  T  E  such  that  (2.1)  is  satisfied.  Moreover, 


l|T||i  +  ||r||o,rc  <  <:^(lbllo.rc  +  IIQIIo  +  \\HN\\o.r.) .  (2.5)1 

The  admissibility  set  Uad  is  defined  by 

Wad  =  {(T,  g)  E  X  L^iTc) :  J{T,  g)  <  oo,  (2.1)  is  satisfied}  .  (2.6) 

Then,  (T,  g)  E  Uad  is  called  an  optimal  solution  if  there  exists  e  >  0  such  that 

J{f,9)<J{T,g)  V  (T,</)eW<.d  satisfying  ||r- fill  +  llff  -  ^llo, r.  <e.  (2.7) 

Based  on  the  previous  lemma,  we  can  show  the  existence  and  uniqueness  of  optimal  solutions. 
Theorem  2.2-  There  exists  a  unique  optimal  solution  (f ,  g)  E  Uad-  ® 


Using  techniques  in  e.g.  [9]  we  may  obtain  the  optimality  condition 


where  $  is  the  solution  of  the  adjoint  state  equation 


Ka(i?,4>)  +  c(u,ii;,4>)-i(/?,r-r<,)p^  =0  (2.9) 

Eliminating  g  from  (2.1)  and  combining  with  (2.9),  we  obtain  the  optimalHy  system 

KaiT,S)  +  c(u,T,S)  +  ji^,S)rc  =  iQ,S)  +  K{HN,S)rs  V  S  €  H},(n)  (2.10) 

and 

Ka(R,^)  +  c(u,R,^)-^(R,T-Td)r^  =0  VRsHhCQ).  (2.11) 

Thus,  the  optimal  state,  i.e.,  the  temperature  distribution  T,  can  be  found  by  solving  the  coupled  system 
(2.10)-(2.11).  which  also  provides  the  optimal  co-state  The  optimal  control  g  can  then  be  deduced  from 
(2.8). 


3.  Finite  Element  Approximation,  Error  Estimates,  and  Iterative  Methods 

In  the  usual  manner,  one  may  construct  finite  element  subspaces  IT’'*  C  Hp{Q)  D  C(n)  parametrized 
by  a  parameter  h  that  tends  to  zero.  (In  practice,  h  is,  of  course,  related  to  a  grid  size.)  We  assume  the 
approximation  property  (see  [3]):  there  exist  an  integer  k  and  a  constant  C  such  that 


inf  115-5'*||i  <C/i"*||51U+i 


V5  6  and  0  <  m  <  ib. 


(3.1) 


A  finite  element  algorithm  for  determining  approximations  of  the  solution  of  the  optimality  system  (2.9)- 
(2.10)  is  then  defined  as  follows:  seek  T*  €  W*  and  6  such  that 

Ka{T\S^)  -h  c(u,r\5'‘)  -h  |($^5'•)^c  ^3  2) 

=  (Q,5'*)-l-/c(f/;v,5'*)r^  V5'‘eW'* 

and 

«a(i^^$'‘)  +  c(u,i^^$'*)-i(i^^^'‘-^d)p^  =0  ^Rl'ew".  (3.3) 

Although  the  optimality  system  is  linear,  the  coupling  of  u  and  $  in  the  two  equations  make  the 
derivation  of  error  estimates  nontrivial.  It  turns  out  to  be  convenient  to  apply  the  Brezzi-Rappaz-Raviart 
theory  (see  [2],  [5],  and  [7])  to  obtain  error  estimates. 


Theorem  3.1-  Lei{T,^)  and{T^ be  the  solutions  of  (2.9)-(2.10),  and  (3.2)-(3.3),  respectively.  Assume 
that  T,  $  €  ^f”*"''^(f2)  n  for  some  I  <m<k;  also  assume  that  (3.1)  holds.  Then, 

<  Cmax  |i,  l|  /l"*(||(3||m-l  +  ll-H'wllrAr.m-l/Z  +  ||7d||rc.m-l/2)  , 

where  C  is  independent  ofh,  6,  j,  T,  and  I 

A  simple  iterative  algorithm  for  solving  (2.10)-(2.11)  can  be  defined  as  follows: 
choose  ; 


for  n  =  1, 2, . . . ,  solve  for  7^”^  from 


/ca(T("),5)+c(u,r(">,5) 


=  +iQ,S)  +  K{Hi,,S)rs  V5  € 


then  solve  for  from 


«  aiR,  $("))  +  c(u,  R,  $(-*))  =  ^{R,  2^”)  -  T^)rc  V  iZ  €  . 


(3.4) 


(3.5) 


Of  course,  ultimately,  this  algorithm  has  to  be  carried  out  in  a  discretized  version,  such  as  one  using  a  finite 
element  method. 

The  convergence  of  this  algorithm  can  be  proved  as  a  result  of  the  observation  that  it  is  effectively  a 
gradient  method  for  the  following  minimization  problem:  find  g  6  I*(rc)  such  that  JC(g)  :=  J{T{g),g)  is 
minimized  where  T{g)  €  ff})(Ci)  is  defined  as  the  solution  of  (2.1). 

Theorem  5.2-  Let  <&("!)  be  the  solution  o/(3.4)-(3.5)  and  {T,i)  the  solution  o/(2.9)-(2.10).  Then, 
T  in  H^o(n)  and  —  $  in  os  n  -»  oo.  I 


4.  Computational  Examples 

Let  Q  C  be  the  unit  square  (0, 1)  x  (0, 1).  Let  F  =  Fc  U  Fp  U  U  To  be  shown  as  in  Figure  1. 


Figure  1.  Computational  domain 


The  finite  element  spaces  are  chosen  to  be  piecewise  linear  elements  on  a  triangular  mesh.  All  the 
numerical  results  make  use  of  the  following  parameters  and  data: 

parameters:  x=  1/0.73;  Re  =  l;  We  =  l;  w  =  1/2 ; 


boundarv  data: 


«yT* 

T=1  on  Fp  ,  "a"  ~  ^  onFjvUFo; 

an 


heat  source: 


5/c  /  3 

Q  =  — Y  (  9’’’*  cos(3xi)  cos'(xy)  —  4x'  sin'(-xr)  cos(2xy) 


+  (1  -  y*)-^  sin(3xz)  cos-(xy) ; 
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velocity  profile: 


u  =  (l-jr,0), 


elastic  extra  stress  profile: 


pressure  profile: 


p  =  4(l-ir)-i. 


In  functional  (1.10),  we  choose 


Td  =  3.5. 


For  the  data  given  above,  the  exact  solution  of  the  uncontroled  problem,  i.e.,  for 


dn 


=  0 


on  Fc, 


is  given  by  T  =  5  sin^(|Ti)  cos^(iry)  +  1. 

We  compare  the  temperature  distribution  in  the  uncontroled  case  with  the  optimal  temperature  distri¬ 
bution  in  the  controled  case  for  which 


dn 


-g  on  Fc , 


where  g  is  the  control  such  that  (1.10)  is  minimized.  Approximations  to  the  optimal  state  cind  co-state 
are  computed  from  (3.2)-(3.3):  the  approximate  optimal  control  p*  is  then  obtsuned  from  (2.8),  i.e.,  g'^  = 
— (l/5)<&*|rc-  All  of  the  computational  results  shown  below  were  obtained  with  the  use  of  a  mesh  size 
/»  =  Of  course,  calculations  with  varying  mesh  sizes  were  performed.  Since  these  merely  verified  the 
error  estimates,  we  do  not  report  on  them  here. 

Specifically,  Figures  2-4  deal  with  the  following  cases: 

1.  Exact,  uncontroled  temperature  and;  Optimal  temperature; 

2.  Exact,  uncontroled  temperature  and;  Optimal  temperature; 

2.  Optimal  boundary  control. 


Figure  2.  Temperature  contours:  uncontroled  and  controled. 
Optimal  boundary  control  on  Fc. 

(Fc  is  the  top  boundary  segment.) 

(Fo  is  the  right  boundary  segment.) 


Figure  3.  Temperature  surfaces;  uncontroled  and  controled. 
(Fc  is  the  top-right  boundary  segment.) 

(Fo  is  the  lower-right  boundary  segment.) 


-6.25 


-9.38 


-12.51 


-15.63 


-18.76  A - - - - - - - - - - - \ 

0.00  0.17  0.33  0.50  0.67  0.83  1.00 

Figure  4.  Optimal  boundary  control  on  Fc- 
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ANALYSIS  AND  FINITE  ELEMENT  APPROXIMATION 
OF  AN  OPTIMAL  CONTROL  PROBLEM  IN 
ELECTROCHEMISTRY  WITH  CURRENT  DENSITY  CONTROLS 


L.  Steven  Houf  and  James  C.  TurnerJ 


Abstract.  An  optimal  control  problem  for  impressed  cathodic  systems  in  electro¬ 
chemistry  is  studied.  The  control  in  this  problem  is  the  current  density  on  the  2mode. 
A  matching  objective  functional  is  considered;  (many  other  objective  functionals  Cein 
be  similarly  treated.)  The  existence  of  an  optimal  solution  is  proved.  The  use  of 
Lagrange  multiplier  rules  is  justified  ^lnd  an  optimality  system  of  equations  is  estab¬ 
lished.  Finally,  a  finite  element  algorithm  is  defined  and  optimal  error  estimates  are 
derived. 


Key  words.  Optimal  control,  impressed  cathodic  system,  electrochemistry,  nonhnear  boundary 
condition,  finite  element  method,  error  estimate 

AMS(MOS)  subject  classifications.  49J20,  49K20,  65K10,  65N15,  65N30 


1.  Introduction 

We  consider  an  optimal  control  problem  for  impressed  cathodic  systems.  A 
typical  example  of  an  impressed  cathodic  system  is  a  metal  container  filled  with 
an  electrolyte.  The  painted  portion  of  the  container  surface  is  usually  treated 
as  instilated.  The  unpainted  part  is  divided  into  cathode  and  anode  which  are 
connected  to  the  negative  and  positive  poles  of  an  electrical  source,  respectively.  By 
adjusting  the  current  density  on  the  anode  we  could  effectively  alter  the  potential 
distribution  on  the  entire  bounding  surface  or  in  the  entire  fiow  domain.  The 
potential  distribution,  of  course,  has  a  direct  effect  on  the  chemical  reaction  process 
occurring  inside  the  fiow  domain,  which  in  turn  has  an  effect  on  the  rate  of  corrosion 
of  the  metal  container.  Thus  the  current  density  on  the  anode  can  be  used  as  a 
practical  control  variable  for  generating  a  desired  potential  field.  This  idea  can  be 
conveniently  formulated  as  optimal  control  problems  for  the  potential  equation  with 
appropriate  boundary  conditions.  Optimal  control  problems  of  this  sort  has  been 
studied  in  [20]  and  [21]  where  the  goal  was  to  match  a  desired  potential  distribution 

tDepartment  of  Mathematics  ajid  Statistics,  Simon  Fraser  University,  Burnaby,  B.C.  V5A  1S6, 
Canada.  Work  supported  in  part  by  the  Natural  Science  and  Engineering  Research  Council  of 
Canada  under  grant  number  OGP-0137436  and  by  a  Simon  Fraser  University  President’s  Research 
grant.  e-mail:  lhou®cs.sfu.ca 

^Department  of  Mathematics,  The  Ohio  State  University,  Columbus,  Ohio  43210-1174,  USA. 
Work  supported  in  part  by  US  Army  Research  Office  under  grant  number  DAALb3-91-G0237. 
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on  the  cathode.  The  models  analyzed  in  [20]  and  [21]  are  essentially  linear.  [14] 
discussed,  mainly  from  an  algorithmic  point  of  view,  several  control  mechanisms 
including  adjusting  the  positions  of  anodes  and/or  the  current  density  on  the  anodes 
in  order  to  best  match  a  desired  potential  on  the  structure  surface;  nonhnear  models 
were  employed  as  well  as  linear  ones.  [15]  analyzed  a  “location  control”  problem, 
z.e.,  the  control  variable  is  the  location  of  anodes,  wherein  nonlinear  models  with 
boundary  conditions  of  polynomial  or  mixed  polynomial-exponential  growth  type 
were  considered.  In  this  article,  we  will  attempt  to  mathematically  analyze  optimal 
control  problems  with  current  density  controls.  The  nonlinear  model  used  involves 
an  exponentially  growing  boundary  condition. 

We  assume  the  electrolyte  occupies  a  physical  domain  G  R  ^  with  a  boundary 
r.  The  domain  is  assumed  to  be  finite  in  this  paper,  although  infinite  domain 
problems  can  be  handled  if  appropriate  decay  rate  at  infinity  is  assumed.  If  C  R^, 
we  will  need  to  work  with  a  non-Hilbert  space  with,  e.^.,  p  =  3;  similar 

results  can  still  be  obtained. 

The  electrical  potential  ^  in  is  governed  by  the  differential  equation 

—  div(<T  grad  ^)  =  0  in  , 

where  the  conductivity  cr  is  a  continuous  function  with  a  positive  lower  bound. 

The  boundary  F  is  divided  into  three  components:  the  anode  F^,  the  cathode 
Fc  and  the  insulated  part  Fq.  On  the  cathode  Vc-,  <f>  satisfies  the  relation 

a^  =  -f(<f>)  onFc, 

where  /  is  an  empirical  function  that  depends  on  the  electrode  materials  (see  [4]). 
In  particiilar  we  will  assume  /  is  given  by  the  Butler- Volmer  function: 

f(<f>)  =  C3[e^‘^  -  (1.1) 

where  C\,  C2  and  C3  are  positive  constants  (see  [4]).  Throughout  this  paper,  /  will 
be  assumed  to  be  defined  by  (1.1).  For  notational  convenience,  we  will  mainly  use 
/(^)  rather  than  the  exphcit  expression  given  in  (1.1). 


FIG.  1.  A  typical  impressed  cathodic  system;  an  electrolyte  container 


connected  with  an  electrical  current  source 


On  the  anode  F^,  we  have  the  boundary  condition 


d(i>  P 

a—  =u  on  Fa  , 
on 

which  corresponds  to  the  specification  of  the  current  density  on  the  anode.  Adjust¬ 
ing  the  current  density  on  F^  amounts  to  treating  u  as  a  control  variable.  On  the 
insulated  part  Fo, 

cr—  =  0  on  Fo  . 
on 

We  are  concerned  with  the  following  optimal  control  problem:  seek  a  state  (f> 
and  a  control  u  such  that  the  functional 


J(M  =  ^  f  (4-M'‘dn  +  ^  f  u^dr, 
^^0  jQ  ^  Jr  A 


(1.2) 


is  minimized  subject  to  the  the  constraint  equations 


—  div(<T  grad  =  0  in  fi 


(1.3) 


on  F^ , 


on  Fo 


and 


on  Fc . 


(1.4) 

(1.5) 

(1.6) 


In  (1.2)  Oo  is  a  desired  potential  distribution  in  Q,  and  eo  and  So  are  positive  con¬ 
stants. 

We  will  mahe  use  of  an  equivalent  variational  formulation  (1.7)  of  the  nonlinear 
boundary'  value  problem  (1.3)-(1.6).  We  will  utilize  Sobolev  spaces  ff^(rA), 

■ff*(Fc)-  ■ff*(Fo)  and  jy*(F).  The  corresponding  norms  on  these  spaces  will  be 
denoted  by,  e.ff.,  ||  •  ||,n,  ||  •  For  details,  see  [1]  and  [8].  A  weah  formulation 

of  (1.3)-(1.6)  is  given  as  follows:  seek  a  G  B'^(Q)  such  that 


/  cr  grad  ^  •  grad  V’ dfi -f-  f  f((f))tpdr=  j  utpdT, 
Jn  JTc  Jta 


VV’  e  H\n). 


(1.7) 


Equation  (1.7)  can  be  formally  shown  to  be  equivalent  to  the  nonlinear  bound¬ 
ary  value  problem  (1.3)-(1.6).  We  mention  that  second  order  elliptic  differential 
equations  with  exponentially  growing  coefficients  were  studied  in,  among  others, 
[9],  [12].  [13]  and  [16].  An  elliptic  equation  with  mixed  Dirichlet-Neumann  type 
boundary'  conditions  that  have  an  exponentially  growing  coefficient  in  the  boundary 
condition  was  studied  [7].  Some  of  the  techniques  in  these  articles  are  useful  for  the 
mathematical  and  numerical  analysis  of  the  state  equation  (1.7). 

We  restate  the  minimization  problem  as  follows: 


seek  a  state  (j)  G  H^{Q)  and  a  u  ^  U  such  that  the 
functional  (1.2)  is  minimized  subject  to  (1.7), 


(1.8) 


where 


U  is  a.  non-empty,  closed,  convex  subset  of  L^(Ta)-  (1-9) 

Now  we  state  a  few  useful  facts.  We  set  a  =  niin{C'i, C'2}.  Then  by  Mean 
Value  Theorem  we  obtain 

/(<^)^  =  Cz(Ci  +  C'2)  >  2a<f>^  V  <l>  (1.10) 

where  o  is  between  {Ci<f>)  and  (— C2^).  It  is  also  easy  to  see  that 

f'{<t>)>2a  W<f>.  (1.11) 


The  norm  on  defined  by 

\\\<f,\\\^=:  cr\gT&d(f>\'^dn  +  2a  j  V<?i€ff^(fi) 

is  equi\'alent  to  the  usual  ^r^(fI)-norm  ||  •  ||i,  i.e., there  exist  constants  p  >  0  and 
7  >  0  such  that 

pWoWI^  f  <T\gTa.d<l>\‘^d{l  +  2a  f  dT  >  7  ||<^||^  (f>  e  H\Q) .  (1.12) 

Jii  Jrc 

A  proof  of  (1.12)  can  found  in,  e.g.,  [17]. 

The  rest  of  the  paper  is  organized  as  follows.  In  §2  we  prove  the  existence  and 
uniqueness  of  solutions  to  (1.7)  so  that  the  constraint  equation  is  well-posed.  In 
§3  we  show  the  existence  of  an  optimal  pair  (<^,  u)  that  minimizes  (1.1)  subject  to 
(1.7).  In  §4  we  justify  the  use  of  Lagrange  multiplier  rules  and  derive  an  optimal¬ 
ity  condition.  In  §5  we  discuss  the  regularity  of  optimal  solutions.  Finally  in  §6 
we  define  a  finite  element  algorithm  for  solving  the  optimality  system  and  derive 
optimal  error  estimates. 

2.  Existence  and  Uniqueness  of 
Solutions  to  the  Constraint  Equations 

We  first  examine  the  existence  of  a  solution  to  the  nonlinear  Neumann  type 
boundary  value  problems  (1.7). 

Lemma  2.1  X  be  a  finite- dimensional  Hilbert  space  whose  scalar  product  is  denoted 
by  (•,•)  and  the  corresponding  norm  |  •  |.  Let  F  be  a  continuous  mapping  form 
X  into  X  with  the  following  property:  there  exists  an  r  >  0  such  that 

iF{4>),(f>)>0  e  X  with  \4>\  =  r  . 

Then  there  exists  an  <f>  &  X  such  that 

F  ((/))  =  0  and 


Proof.  See  [Raviart],  pp.  279.  I 


Lemma  2.2  Assume  <f>  6  and  s  >  0.  Then  €  L^(r).  Moreover,  there 

exists  a  constant  k,  independent  of  (j>,  such  that 


J  rfr  <  1  +  |r|  +  |r| 


<  oo , 


where  |r!  is  the  measure  ofT. 


Proof.  Let  (j)  €  and  5  >  0  be  given.  A  Sobolev  embedding  theorem 

implies  o  G  Using  embedding  results  for  Orlicz- Sobolev  spaces  (see  [1], 

[10]  and  [18])  (recall  Q,  C  R^),  we  have  La^T)  where  the  iV-function 

A(t)  =  —  1.  Thus  there  exists  a  constant  k  >  0  such  that 


PIUA(r)  = 


r  /-Ml! 

inf  <  A:  :  y  ^ 


dr  <  1  >  <  /c||<^||i  <  oo. 


Hence  for  each  sufficiently  e  >  0,  the  constant  k  =  ||<^||LA(r)  +  ^  satisfies 

J  [e^  -  1]  dr  <  1  so  that  j  e^  dT  <  1  +  |r| . 

We  set  M  =  sk^.  By  an  elementary  calculation  we  can  show  that 

e«<e^  V|r|>M. 


We  set  K  =  =  e"  *  <  oo.  Then 

f  e"l^ldr=  /  e"l^ldr+  /  e'l-^'dr 

Jt  J{x€r-.\,t>(x)\>M]  y{x€r:|^(x)|<Af} 

<  j  <ir  +  JiT  |r|  <  1  +  |r|  +  |p| 

<  1  +  in  +  |r|<tx.. 

Letting  e  — >  0  yield  the  desired  result.  I 


Lemma  2.3  Assume  {<^n}  C  Li^iTc)  is  a  sequence  such  that  (j>n  (f>  a.e.  on  Tc 
and 

r 

/„  dr  <  J5  V  n 

c 


f  f(4>n)<i>n<Ii^  <  B  Vn  (2.1) 

Jrc 

where  f  is  defined  by  (1.1)  and  B  >  0  is  a  constant  independent  of  n.  Then 

[  /(^)<^dr  <liminf  /  f{<l>n)<i>ndr 
JTc  JTc 

and 

hm  /  |/(?in)-/(?i)|dr  =  0. 
oo  Jvn 
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Proof.  The  proof  follows  the  ideas  of  [13],  pp.  21-22.  Since  /  is  continuous  and 
(j)n  o  a.e.  on  Tc,  we  deduce  that  f(<f>n)  —*  a.e.  on  Fc-  Note  that  f{(j>)<f>  >  0 
on  Fc  so  that  we  may  use  Fatou’s  Lemma  to  obtain 

f  f(4>)<f)dT  <liminf  f  dF  <B. 

Jtc  JTc 

Hence  fi<f>)4>  €  By  setting  K  =  sup|jj<i  |/(a:)|  we  easily  conclude  from  the 

identity 

\m\  =  \t\-^mt  vt/o ,  (2.2) 


that 

Thus 


\m\<mt+K  vieR. 


\m\<\m\\<f>\+K  onFc, 


i.e.,  f{o)  €  L^(Fc).  Utilizing  (2.2)  again,  we  deduce  that  for  each  ^  >  0  and  for 
a.e.  X  €  Fc,  we  have  either 


\(t>n\<^  ^  or  \f{<f>)\<Sf{4>n)(i>n 


SO  that 


\f{<t>)\  <C6+  Sf(<f>n)<f>n  on  Fc  , 


where  Cs  =  sup  |/(a:)|.  For  every  measurable  subset  5  C  Fc  we  have 


/  \f{4>n)\dr<Cs\S\  +  6  f  f(<l>n)4>ndr. 
JTc  JTc 


Equation  (2.1)  implies 


2B 


J  fi4>n)<l>ndT  < 
for  n  greater  than  some  No  >  0.  Thus 

J^\f((l>n)\dr<C6\S\  +  2B6  Vn>iVo, 

where  15|  is  the  measure  of  S.  Hence,  the  sequence  of  functions  {f{(l>n)}  has  equi- 
absolutely  continuous  integrals.  By  Vitali’s  Convergence  Theorem, 

Km  f  \f(4,„)-fW\dr  =  0.  I 


Theorem  2.4  Assume  u  €  L^{Va)-  Then  there  exists  a  unique  4>  6  Ff^(fi)  that 
satisfies  (1.7).  Furthermore,  (f>  satisfies  the  estimate 


where  /3  and  7  are  constants  independent  of  <f>. 
Proof.  We  first  establish  an  inequality 


/  a  grad  <f>  •  grad  (f>dQ,  +  j  f{(f>)<f>dr  —  j  ucjxir  > 

Jq  Jtc  Jt'a 

>0  V  E  H^{Q,)  with  ||<^||i  =  r 


for  some  r  >  0.  Using  (1-10)  we  have  that 


<fr  >  O' 


<ir 


w<f>e  H\Q). 


Hence  it  follows  from  (1.12)  that 


/  cr  grad  ^  •  grad  dfl  +  f  dT  >  j\\(f>\\i 

Jsi  Jfc 

Using  Cauchy- Schwartz  inequality  and  trace  theorems  we  obtain  the  estimate 


where  3  is  a  positive  constant.  By  combining  the  last  two  estimates  we  deduce  that 
for  r  =  ^  {Iloilo, r>t  -I- 1}  >  0  we  have 


/  <T  grad  ^  •  grad  ^  -|-  /  f{<j>)<f>dT—  j  ucpdF  > 

Jii  JTc  Jta  (2.4) 

>0  'i  (j>  E  H^{Q)  with  ||(^|(i  =  r . 

Since  is  separable,  we  may  choose  a  countable  basis  of  We 

set  Xn  =  span{V’i,-  •  •  ,  V’n}-  The  inner  product  and  norm  on  each  Xn  is  defined 
by  that  of  restricted  to  X„.  We  introduce  the  mapping  Fn  :  Xn  —*■  Xn  as 

follows.  For  each  (j>  €  Fn(^)  E  Xn  is  defined  by 

{Fn{<f>).uj)  =  I  <7  grad  •  grad  V’ j  +  /  /(^)i/’j  dT  —  j  uipjdT  l<j<n. 

Jn  Jtc  Jta 

It  follows  from  (2.4)  and  Lemma  2.1  that  the  finite  dimensional  problem 

I  grad  •  grad  V’ dfi -f  /  /(^„)V’dr  =  j  u^dT,  V  V’ €  (2.5) 

dn  Jtc  Jta 

has  a  solution  E  Xn  with  a  bound 

ll<^Ti||i  <  ^  {||«lio,rA  +  1}  •  (2.6) 
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We  can  extract  a  subsequence  of  still  denoted  by  that  converges  weakly 
to  some  6  €  as  n  0.  Then  {<f>n}  also  converges  weakly  in  by 

a  trace  theorem.  Thus  {<i>n}  converges  strongly  in  I'^(r)  by  compact  imbedding. 
This  in  turn  implies  a  subsequence  satisfies  (f>n  <f>  a.e.  on  F.  By  setting  if)  =■  (f>n 
•  in  (2.5)  we  obtain 


/  <  ||u||o,rx  llfl^n||o,rx 

JTc 


<  ^ll“llo,rx  ll<?^n||i  <  c  l|u||o,rA  -  {||“llo,rA  +  1} 


By  Lemma  2.3  we  have  that 


dr  <  liminf 

n— ►oo 


dr 


and 


lim  /  \f{<i>n)  -  f{4>)\ dT  =  0 . 

Jtc 

By  passing  to  the  limit  in  (2.5)  for  each  r/)  6  C°°{Q)  we  see  that 

f  <7  grad  ^  •  grad  V’ dn  +  /  dT  =■  j  ut/^dT,  V  G  , 

Jci  Jr  c  Jr  A 


(2.7) 


#  We  next  prove  that  this  </>  is  a  solution  to  (1.7).  For  each  we  may 

choose  a  sequence  {V’*}  C  C°°(fl)  such  that  |j^  —  V’Jklli  — >•  0  as  — >•  oo.  Using  (2.7) 
we  have 


j  (T grad •  grad V’fc  dfl  +  j  /(^)^jt  dT  =  j  ut/’t  dT ,  Vfc.  (2.8) 

J  Q  JTc  a 


Lemma  2.2  implies  €  L^{Tc)  and  e  €  L^{Tc)  so  that  |/(<^)P  cfT  <  oo. 
Hence. 


V’jfc)dr 


<  ||/(<^)||o,rc  IIV’  -  ^fc||o,rc  <  C'll/(<^)llo,rc  U  -  Mi 


so  that 


/(</.)(^  -  ^fc)dr  0 


as  k 


oo . 


Thus  we  may  pass  to  the  limit  in  (2.7)  to  show  (1.7)  holds.  The  estimate  (2.3) 
follows  easily  from  (2.6). 

To  answer  the  question  of  tmiqueness,  we  assume  ^  and  ^  axe  two  solutions  to 
(1.7).  Then  we  have 


f  (7grad(<^  -  ^)  •  gradV’dO  +  /  [/(<^)  - /(^)]^  dT  =  0 ,  Vt^GJEf^(O). 
Jn  Jrc 
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Setting  V  ■=  (j>  —  (f>,  we  see  that 

/  cr|grad((?!>- f  [f(<j))-f(^)]{<j>-^)dr  =  0. 
Jn  JTc 

Using  Mean  Value  Theorem, 

[  (7|grad((^  -  +  /  f  dT  =  0 

J(1  JTc 


for  some  ^  between  ^  and  <f>.  Using  (1.11),  we  see  that 

/  a\  grad((^  —  (^)p  dO,  +  a  f  ((/>  —  dT  <  0 . 
Jq  JTc 


Hence  we  deduce  that  grad(<^  —  (^)  =  0  in  and  ((^  —  (^)  =  0  on  Fc-  This  in  turn 
implies  (d  —  ^)  =  0  in  fl,  t.c.,  uniqueness  holds.  I 


3.  Existence  of  An  Optimal  Solution 

Having  shown  that  the  constraint  equation  (1.7)  is  well-posed,  we  are  now  pre¬ 
pared  to  study  the  existence  of  an  optimal  solution  (<l>,  u)  that  minimizes  the  func¬ 
tional  (1.2)  subject  to  (1.7).  We  introduce  the  admissible  set 

l^ad  =  {(<^5^)  e  H^(Q.)  X  U  :  (<f>,u)  satisfies  (1.7)}, 

where  U  is  given  by  (1.9). 

Theorem  3.1  There  exists  a  (^,u)  €  x  U  that  minimizes  (1.2)  subject  to 

(1.7). 


Proof.  Theorem  2.4  implies  an  element  {(f>,  u)  E  Had  exists  such  that  v7(^,  u)  <  oo. 
Let  {(<?!»„,«„)}  C  Had  be  a  minimizing  sequence,  z.e.. 


Vila.  J{(f)n,Un)=  inf  J{<t>iu)  (3.1) 

n— oo  (<^,u)eMod 


and 


I  cr grad -grad V’dfl-H  f  f{4n)i’dr=  f  Unil>dr,  VV’Gif^(fi).  (3.2) 
Jii  JTc  Jta 


Using  (1.2)  and  (3.1)  we  deduce  {un}  is  bounded  in  L^(r^).  Then  (2.3)  implies 
{||<^n||i}  is  bounded.  Hence  we  may  extract  a  subsequence  {{(l>n,Un)}  such  that 

<f)n  (f>  in  H  {^)  and  ^  u  in  L  (Ta)  • 


Furthermore,  trace  theorems  implies  <f>n  ^  4>  'va.  L'^{Tc)\  this  in  turn  implies  <j>n  ^  <i> 
a.e.  on  Fc  (after  extracting  subsequences  if  necessary).  By  setting  ^  in  (3.2) 
we  obtain 


i“nl|o,rA  ll<^n||o,r^  <  /?  ||«n||o,rA  • 
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Hence  we  deduce 


/  fi<f>n)<f>n  <  M 

JTc 

where  AI  is  a  constant  independent  of  n.  By  Lemma  2.3, 

/  /((^)(^dr  <  liminf  /  f((f>n)4>ndT 
JTc  Jtc 

and 

lim  /  \fi<f>n)  -  /(<^)|  =  0.. 

By  passing  to  the  limit  in  (3.2)  for  tp  €  C°°(Q,)  we  obtain 

I  a  gradp  •  gr&drp  dQ  +  j  f{^)‘tpdT—  j  uipdF,  V  V’ G  . 

Jq  Jt  c  Jt  a 

Then  using  the  denseness  of  in  H^(Q,)  and  the  fact  that  f{p)  €  Z’^(rc'),  we 

obtain 

j  O’ grad  •  grad  V’ df)  +  /  f($)ipdr=  j  uip  dT ,  V  . 

Jq  Jtc  X 


Thus  (o.u)  G  Uad.  Finally  using  the  weak  lower  semi-continuity  of  J {•,'),  we 
conclude  that  (<^,  u)  is  indeed  an  optimal  solution,  i.e., 

J{^,u)=  inf  J{<i>,u).  I 

{<i>,u)£Uad 

4.  Lagrange  Multiplier  Rules 

In  this  section  we  will  attempt  to  characterize  optimal  solutions  whose  existence 
has  been  established  in  §3.  Since  the  constraint  equation  has  a  unique  solution  for 
each  given  control  u,  the  state  ^  is  a  well-defined  function  of  u.  Since  U  given  in 
(1.9)  is  not  necessarily  an  open  set,  the  mapping  u  ^  is  in  general  not  differen¬ 
tiable.  Although  other  approaches  available,  it  turns  out  to  be  convenient  to  use 
the  Lagrange  multiplier  riile  to  turn  the  constrained  minimization  problem  (1.8) 
into  an  unconstrained  one  and  derive  an  optimality  condition. 


We  first  quote  the  following  abstract  theorem  concerning  the  existence  of  Lagrange 
multipliers  for  minimization  problems  on  Banach  spaces  (see,  e.g.,  [19]): 

Lemma  4.1.  Let  B\  and  B2  be  two  Banach  spaces,  and  U  an  arbitrary  set.  Suppose 
J  is  a  functional  on  B\  x  U ,  K  a  mapping  from  B\  x  U  to  B2.  Assume  (p,u)  € 
B\  xU  is  a  local  minimum  for  the  constrained  minimization  problem: 

min  J'(<j),u)  subject  to  )C(<l>,u)  =  0,  (4-1) 


i.e.,  there  exists  an  open  neighborhood  Q  of  ^  in  B\  such  that 

J (^,  u)  <  J {<i>,  u)  W  {<l),u)  E  Q  X  U  satisfying  }C{<f>,  it)  =  0 . 


Assume  further  the  following  conditions  are  satisfied: 

(A)  for  each  u  £  U,  and  <}>  K{<f>^u)  are  Frechet- differentiable  in 

0: 

(B)  for  every  (f>  ^  Q,  ui,U2  E  U  and  a.  E  [0,1],  there  exists  a  Ua  =  Uo,(<f>,Ui,U2)  E 
U  such  that 

/C(^,Ua)  =  a)C(<l>,ui)  +  (1  —  a))C{4>,U2) 

and 

J{<i>-,Ua)  <  OiJ{<j>,u{)  +  (1  -  OL)J{<i),U2)  ; 

(C)  the  algebraic  sum  of  u)Bi  +  AC(^,  U)  contains  a  neighborhood  of  0. 

Then  there  exists  a  X  E  B2  such  that 


and 


mmC($,u,X)  =  C(^,u,X), 

u£U 


where  C{f),u,X)  =  J'(<j>,u)  —  {X^lC((f>,u))  is  the  Lagrangian  of  the  constrained  min¬ 
imization  problem  (4.1). 

Proof;  See  [19].  I 


We  mil  fit  our  optimization  problem  (1.8)  into  the  above  abstract  framework. 
We  define  the  Banach  spaces  Bi  =  H^(Q)  and  B2  =  Let  U  be  given  by 

(1.9).  The  (generalized)  nonlinear  constraint  1C  :  BixU  — +  B2  is  defined  as  follows: 
AC(^,  u)  =  I  for  (^,  u)  E  B\  X  U  and  f  €  152  if  and  only  if 


(/,V’)rx  =  /  O’ grad  •  grad  V’df)  +  /  f(^)ifdr 
JQ  JTc 


-  [  uxfdT  WxfE  H\n). 

Jta 


(4.2) 


We  easily  see  that  (1.8)  is  equivalent  to 

find  (</»,  u)  E  Bi  xU  such  that 

J{(t>,u)  —  inf  :  (^,u)  E  Bi  x  U,  K{^,u)  =  0}  . 

Let  the  Frechet  derivative  of  J"  and  K  with  respect  to  <f>  be  denoted  by  DijiJ'  and 
respectively.  Let  the  Frechet  derivative  of  J  and  )C  with  respect  to  u  be 
denoted  by  Duff  and  respectively.  D^)C{(j>,u)  E  C{Bi,B2)  is  defined  as 

follows.  u)  •  f>  =  l  iox  E  B\  and  f  €  B2  if  and  only  if 

(r.u)r^=  /  O’ grad  ^  •  grad  ^  dfl  +  /  f'{(l>)^ipdr  VV’€if^(fl).  (4.3) 

Jq  Jtc 


In  order  to  apply  Lemma  4.1,  we  need  to  verify  conditions  (A)-(C).  We  begin  with 
the  verification  of  (C). 


Lemma  4.2.  Assume  ti)  €  Had  optimal  solution  to  (1.7).  Then  the 

operator  Dfj,K,{(f>^u)  is  onto  from  B\  to  B2- 

Proof:  Using  (1.11)-(1.12)  we  easily  obtain  the  coercivity  for  the  bilinear  form 


{o.xj})^  f  cr  grad  ^  •  grad  ^  +  /  f' dT  'i  <i> ,  if  ^  {€l) . 

Jq  Jtc 

Then  Lax-Milgram  Lemma  implies  that  for  any  /  €  there  exists  a  imique 

^  that  solves  (4.3).  I 

Now  we  are  prepared  to  derive  an  optimality  condition. 

Theorem  4.3.  Assume  (^,  u)  €  Had  ^  an  optimal  solution  to  the  minimization 
problem  (1.8).  Then  there  exists  a  X  G  H^(Q)  such  that 

f  <7  grad  V)  •  grad  A  dfi  +  f  f\^)xj)\dr  =  —  f  —  <l>o)'if  dO, 

Jn  Jtc 

(4.4) 

and 

f  {6ou  +  X)(u-u)dr>0  \/ueU.  (4.5) 

JTa 


Proof:  Condition  (A)  in  Lemma  4.1  is  obviously  satisfied  for  the  present  setting. 
The  convexity  in  u  for  the  constraint  is  readily  verified  since  the  control  variable  u 
enters  the  constraint  equation  in  a  linear  manner  and  the  control  set  U  is  convex. 
The  convexity  in  u  for  the  functional  is  also  easily  seen  from  the  convexity  of  the 
mapping  u  i— »■  ||u||2  p^.  Thus  Condition  (B)  is  verified.  The  validity  of  Condition 
(C)  is  established  in  Lemma  4.2.  Hence,  by  Lemma  4.1,  there  exists  a  A  6  = 

that  satisfies 

u),  tP)  -  (A,  u)--tP)=0  V  V’  e  (f2)  (4.6) 


and 


where 


min£(<?i,  u.  A)  =  C{<f>,  u,  A) , 

u£U 


£,{<!>, u.  X)  =  J'{(j>,u)—\  f  <7  grad  •  grad  A  dfi  +  j  f(<l))Xdr  —  j  uAdT 

Usi  JTc  Jta 


(4.7) 


•  (4.8) 


Using  (1.2)  and  (4.3),  which  are  the  defining  equations  for  J{^,u)  and 
respectively,  (4.6)  can  be  rewritten  as 

—  I  i<?  —  <f>o)'>P  dO.  —  f  <7  grad  V’ •  grad  A  dJJ  —  /  f'{^)ipXdr  =  0  V^€B^(ft) 

^0  JQ  Jq  Jrc 


which  is  clearly  equivalent  to  (4.4).  Equations  (4.7)-(4.8)  implies  that,  for  ail  u  G 

^0 

'rc 


-  f  (^-ci>o)^dQ  +  ^  f  u^dT 
^0  Jq  ^  JYc 

—  /  <T  grad  •  grad  A  +  /  f(^)Xdr—  I  uXdT 

Un  Jtc  Jta 

cl/ +  ^  /  ^,2^ 

^0  Jn  2  Jpc 

—  I  <T  grad  ^  •  grad  A  dfi  +  /  /(<^)Adr—  j  uXdT  , 

Un  Jtc  JTa  J 


i.e.. 


j  (u  —  ■u)(u  +  u)  «ir  +  /  (u  — ■u)Adr>0  'iu^U. 

Jtc  JTa 


(4.9) 


b. 

2  ^rc 

Given  einy  w  £U  and  e  €  (0, 1),  we  set  u  =  (1  —  e)u  +  ew  in  (4.9)  (note  u  gU  due 
to  the  convexity  oi  U)  and  obtain 


b  j  e  (w  —  u)(ew  —  €u  +  2u)  dr  +  j  e{w  —  u)Xdr>0  w  £  U  . 

-  he  JTa 

Thus  (4.5)  follows  by  dividing  the  last  inequality  by  e  and  then  letting  e  — O'*’.  I 

5.  An  Optimality  System  And  The  Regularity  Of  Its  Solutions 

In  the  sequel  we  will  treat  the  special  case  U  =  L‘^{Ta)-  From  (4.5)  we  easily 
obtain 

u  =  — -^A.  (5.1) 

do 

From  (4.4),  (5.1)  and  the  original  constraint  equation  (1.7),  we  form  the  following 
system  of  equations  (dispensing  with  the  hat  notations  to  denote  optimal  solutions): 

/  (7 grad ^  •  grad ^ dfl  +  /  /(^)^dr  =  — 1  /  XipdT,  V  xj;  £  H^{Q,)  (5.2) 

Jn  Jtc  Jva 


and 


/ 


grad  A  •  grad  a;  dO,  + 


f  f{<f>)Xu;dr 

JTc 

—  —  / (^  ~  (l>o)u;dQ ,  V  a;  6  H^{Q,) . 

eo  Jn 


(5.3) 


This  system  of  equations  will  be  called  the  optimality  system. 

Integrations  by  parts  may  be  used  to  show  that  the  system  (5.2)-(5.3)  constitutes 
a  weak  formulation  of  the  problem 


div(cr  grad  ^)  =  0  in  fl , 


(5.4) 


dS  1  dS  d<f> 

cr-^  =  Fa  ,  <7—  =  0  on  To  and  cr—  =  on  Fc  ,  (5.5) 


13 


—  div((T  grad  A)  =  — (^  —  <^o)  infi, 
Co 


dx 


a—  =0  on  Ta  U  To 
on 


and  =  -f{(l>)X  on  Tc  , 


(5.6) 

(5.7) 


Now  we  examine  the  regularity  of  solutions  of  the  optimality  system  (5.2)-(5.3), 
or  equi^•alently,  (5.4)-(5.7). 


Theorem  5.1  Suppose  that  Q,  C  is  convex  or  of  class  Assume  (^,  A)  € 

H^{Q)  X  fr^(fi)  is  a  solution  to  the  optimality  system  (5.2)-(5.3),  or  equivalently, 
•  (5.4)-(5.7),  then  we  have  that  (<f>,X)  G  x  for  r  G  [l,oo). 

Proof:  Since  A  G  Lemma  2.3  implies  /(<^)  G  I-’'(rc)  and  f'{(t>)  G  L''{Tc) 

for  all  r  G  [l,oo).  We  infer  from  trace  theorems  that  A  G  L'^iVc)  for  all  9  >  1. 
Hence  we  have  G  L^(T')  and  ^  jL®(r)  for  each  q  >  1.  By  applying 
^  elliptic  regularity  results  to  equations  (5.4)-(5.7),  we  obtain  (f  G  W^/^’*(f2)  and 

A  G  W^/^’®(f2)  for  each  q  >  1. 


Remark  In  general  the  possible  discontinuity  of  the  normal  derivative  on  the  inter¬ 
section  ofVc,  To  and  prohibits  us  from  obtaining  further  regularity.  However, 
if  (f>  and  A  vanish  on  the  entire  intersection  of  Fc,  Fq  and  F^,  then  we  could  in 
fact  show  that  (j>  G  C^(f2)  D  C7(f2)  and  A  G  ^7^(12)  Pi  C(f2),  i.e.,  <f>  and  A  are  in  fact 
classical  solutions  of  the  optimality  system.  Also,  H^{Q) -regularity  for  (f>  and  A  is 
expected. 


^  6.  Finite  Element  Approximations 

6.1  Finite  Element  Discretizations.  A  finite  element  discretization  of  the  opti¬ 
mality  system  (5.2)-(5.3)  is  defined  in  the  usual  manner.  For  simplicity  we  assume 
the  domain  is  a  convex  polygon.  We  first  choose  families  of  finite  dimensional  sub¬ 
spaces  1  *  C  H^{^)  satisfying  the  approximation  property:  there  exists  a  constant 
^  C  and  an  integer  k  such  that 


^  v  ^  ,  I  <  m  <  k .  (6.1) 

One  may  consult,  e.g.,  [3]  or  [6]  for  a  catalogue  of  finite  element  spaces  satisfying 
(6.1).  Then,  we  may  formulate  the  approximate  problem  for  the  optimality  system 
#  (5.2)-(5.3):  seek  (j>^  G  and  A*  G  such  that 


/ 


grad  •  grad  V’*  dO,  -f 


/  f{<f>^)^^dT 

JTc 

=  X^^^dT,  VV’*  € 

oo  Jr. 


(6.2) 


and 


I  crgrad  A^  •  grado;*  dfi  4-  [  /'(^*)A*u;^  dF 
Jn  Jtc 

=  —  /(<^*-(^o)a;*dfi,  Vu;*gF' 
^0  Jn 


(6.3) 
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6.2.  Quotation  of  Brezzi-Rappaz-Raviart  Approximation  Theory. 

The  error  estimate  to  be  derived  in  Section  6.3  malces  use  of  results  developed  by 
Brezzi,  Rappaz  and  Raviaxt  (see  [5],  also  [7]  and  [8])  concerning  the  approximation 
of  a  class  of  nonlinear  problems.  Here,  for  the  sake  of  completeness,  we  will  state 
the  rele^•ant  results,  specialized  to  our  needs. 

The  nonlinear  problems  considered  in  [5]  (also  [7]  and  [8])  are  of  the  type 

F(^)  =  +  TG(^)  =  0  (6.4) 

where  A”  and  Y  are  Banach  spaces  and  T  €  JC.(Y ;  A),  G  is  a  G^  mapping  from  X 
into  Y.  A  solution  xjj  to  the  equation  i^(V’)  =  0  is  called  a  nonsingular  solution  if 
we  have  that  F^tp)  is  an  isomorphism  from  X  into  X.  (Here,  F'(-)  denotes  the 
Frechet  derivative  of  F'(-).) 

Approximations  are  defined  by  introducing  a  family  of  finite  dimensional  sub¬ 
spaces  A'*  C  X  and  for  each  >  0  an  approximating  operator  €  C{Y;X^). 
Then,  we  seek  €  X^  such  that 

F^(xP^)  =  rp^ +T^G{ip^)  =  0.  (6.5) 

We  will  assume  that  there  exists  another  Banach  space  Z,  contained  in  F,  with 
continuous  imbedding,  such  that 

G'(V>)  €  £(A; Z)  ^i>eX. 

Concerning  the  operator  T*,  we  assume  the  approximation  properties 

lim||(r'‘-r)y|U  =0  Vy€F 

ft— ►O 

and 

lim||(T*-T)|U(Z;X)=0.  (6.8) 

ft— ^0 

Note  that  (6.6)  and  (6.8)  imply  that  the  operator  G'(V’)  G  £(A;A)  is  compact. 
Morevover,  (6.8)  follows  from  (6.7)  whenever  the  imbedding  Z  C  F  is  compact. 

We  ceui  now  state  the  first  result  that  will  be  used  in  the  sequel.  In  the  statement 
of  the  theorem,  G"  represents  the  second  order  Frechet  derivative  of  G. 

Theorem  6.1  Let  X  and  Y  be  Banach  spaces.  Assume  that  G  is  a  second  order 
Frechet  differentiable  mapping  from  X  into  Y  and  that  G"  is  bounded  on  all  bounded 
sets  of  X.  Assume  that  (6.6)-(6.8)  hold  and  that  ip  is  a  nonsingular  solution  of 
(6.4).  Then,  there  exists  a  S  >  0  and  an  ho  >  0  such  that  for  h  <  ho,  there  exists  a 
unique  €  X^  satisfying  ip^  is  a  nonsingular  solution  of  {6.5)  and  ^  h. 

Moreover,  there  exists  a  constant  G  >  0,  independent  of  h,  such  that 

||V’*->/-|U  <C||(T‘-T)G(«/>)|lx  I 

For  the  second  result,  we  need  to  introduce  two  other  Banach  spaces  H 
such  that  W  C  X  C  H,  with  continuous  imbeddings,  and  assume  that 

for  all  w  €  W,  the  operator  G'(w)  may  be  extended  as  a  linear 
operator  of  C{H',Y),  the  mapping  w  — >  G'{w)  being  continuous 
from  W  onto  C{H-,Y). 


(6.9) 

andW, 

(6.10) 


(6.6) 

(6.7) 


We  also  suppose  that 


(6.11) 


Ji™  -  T\\c(Y-,H)  =  0  ■ 

A — ►O 

Then  we  may  state  the  following  additional  result. 

Theorem  6.2  Assume  that  the  hypotheses  of  Theorem  6.1  hold  and  that  (6.10)  and 
(6.11)  hold.  Assume  further  that 

is  an  isomorphism  of  H  .  (6.12) 

Then,  for  h  <  h\  sufficiently  small,  there  exists  a  constant  C ,  independent  of  h, 
such  that 

||i>'‘-V>||H<c||(r‘-r)G(0)||H  +  ||v.‘-v-|l5c.  I  (6.13) 

6.3  Error  Estimates  for  the  Approximations  of  Solutions  of  the  Optimal¬ 
ity  System. 

In  order  to  derive  error  estimates,  we  begin  by  recasting  the  optimality  system 
(5.2)-(5.3)  and  its  discretization  (6.2)-(6.3)  into  a  form  that  fits  into  the  framework 
of  Brezzi-Rappaz-Raviart  theory  summarized  in  §6.2. 

We  define 

X  =:  X  , 

Y  =  X  X  , 

Z  =  I2(r)  X  X  L'^iV) 

and 

X*  =  y*  X 

where  denotes  the  dual  space  of  Note  that  using  Sobolev  imbedding 

theorems,  Z  CY  with  a  compact  imbedding. 

Let  the  operator  T  G  C{Y-,X)  be  defined  in  the  following  manner:  T{(^,T),d)  = 
{(f),  A)  for  (C,  T],d)  £Y  and  {(f),  A)  G  A  if  and  only  if 


I  cr  grad  •  grad  V’ dfl  +  o  /  (ftp  dT  =  {(,if))r  ,  VV>Gff^(fi)  (6-14) 

J  n  Jtc 


and 


j  (Tgrad  A  •  grades  dfi  +  a  j  Xu)  dT  =  {i],u)  +  {9,u)t  ,  Vu)^H^{Q).  (6.15) 

o/i2  Jt  c 


Clearly,  (6.14)-(6.15)  consists  of  two  uncoupled  elliptic  equations  with  mixed  Robin- 
Neumann  type  botmdary  conditions  and  T  is  its  solution  operator. 

Analogously,  the  operator  G  C{Y;X^)  is  defined  as  follows:  T^{(,t],6)  = 
(^^,  A*)  for  {(^,r),6)  G  Y  and  (<^^,A*)  G  X^  if  and  only  if 

/  (T  grad (/>*  •  grad dfl  +  a  /  (p^tp^  dT  =  {(^,if)^)r  ,  V  G  V*  (6.16) 

Jn  Jtc 


and 


/<Tg^adA'‘•gradu;^d^^+a  /  dF  =  (r/,a;'*)+(0,a;*)r  ,  Vw'*  €  (6.17) 

V  •'Pc 

Clearly,  (  6.16)-(6.17)  consists  of  two  discrete  Poisson-type  equations  that  are  dis¬ 
cretizations  of  the  equations  (6.14)-(6.15);  also,  is  the  solution  operator  for 
these  two  discrete  equations. 

By  the  well-known  results  concerning  the  approximation  of  elliptic  equations 
(see,  e.g..  [3]  or  [6]),  we  obtain: 

||(T-T'‘)(C,r?,^)||x^0  ash-^O,  (6.18) 

for  all  {Q.rj^d)  G  Y  and,  in  addition,  if  G  x  then 

||{T-T'')(C,./,«)||x  <CA”||T«,,,«)||H™+.(n,xH™+.(!!).  (6.19) 


Also,  because  Z  CY  with  a  compact  imbedding,  we  have  that 


ll(T-T'‘)|U(Z;X)^0  asA-.0. 

(6.20) 

Next,  we  define  the  nonlinear  mapping  G  :  X  — »•  F  as  follows:  G((^,  A) 
for  {(f),  A)  €  X  and  (C,  Tj,6)  ^Y  if  and  only  if 

II 

II 

-t- 

t 

f  (f{<j>)  -  a(f>)iT  dT  Vtt  G  iy^/^(r), 

he 

(6.21) 

Co  J 

f  {<f>  -  <j>o)  u dQ.  Va;G^ri(fi) 

(6.22) 

and 

A 

{0,r)T=  if  (<!>)- a)XT  dr  VrGfr'/"(r). 

Jrc 

(6.23) 

(6.21)-(6.23)  is  equivalent  to 

c=  < 

onP^i; 

f{4>)-a<p  on  Pc; 

,  0  on  To  , 

(6.24) 

77  = 

-—{(!>  —  <f)o)  in  fi 

Co 

(6.25) 

and 

^  r(m-«)A  onrc. 

(0  onPoUr^i. 

Recall  /(o)  =  Cz(e^^'^  —  so  that  f'(<f>)  =  C'3(Cie^i'^-l-C2e^^'^).  Using  Lemma 

2.2  and  trace  theorems  we  infer  that  if  (^,  A)  G  x  ff^(fl),  then  for  all  5  >  1, 


17 


(P\t  €  L'^{T),  A|r  €  ■£'’(r),  f{<^)  €  -£'®(r)  and  f'(4>)  G  i®(r).  Hence  we  see  that  the 
triplet  (C,  ?7,^)  defined  by  (6.24)-(6.26)  is  indeed  in  Y,  z.e.,  G  is  well-defined. 

It  is  easily  seen  that  the  optimality  system  (5.2)-(5.3)  is  equivalent  to 

(^,X)  +  TGi(f>,X)  =  0  (6.27) 

and  that  the  discrete  optimality  system  (6.2)-(6.3)  is  equivalent  to 

(4>^,X^)  +  T^G{<f>^,X^)=0.  (6.28) 

We  have  thus  recast  otur  continuous  and  discrete  optimality  problems  into  a  form 
that  enables  us  to  apply  the  theories  of  §6.2.  It  remains  to  verify  the  hypotheses  in 
Theorem  6.1.  This  will  be  the  task  of  the  next  two  propositions. 


Proposition  6.3  The  operator  G  :  X  Y  defined  by  (6.21)-(6.23)  is  second  order 
Frechet  differentiable.  Furthermore,  (6.6)  holds  and  G"  is  bounded  on  all  bounded 
sets  of  X. 

Proof.  In  showing  the  differentiability  of  G,  the  linear  terms  appearing  in  the 
definition  of  G  does  not  pose  any  difficulty.  Furthermore  the  nonlinear  terms  in 
(6.21)  and  (6.23)  can  be  dealt  with  in  a  similar  way.  For  clarity,  we  will  only 
analyse  the  differentiability  of  the  nonlinear  term  t  f'((f))XT  dT.  We  define 

a  mapping  Q  :  X  by  {Q(f>,X),T)  =  /'(^)Ardr  for  all  (<^,  A)  G  X 

and  r  6  For  each  given  ((f),  X)  £  X  we  have  that 


(6.29) 


(^Q(<f,  +  S(f,,X  +  6X)-Q(<l),X),T) 

-f  f"((f>)(S(f>)XTdr- [  f'((t>)(sx)Tdr 
Jtc  Jrc 

=  f  [f'((f, + 8(f>)  -  n<j>)  -  r\<t>m)]  xr  dr 

JTc 

+  /  [f'{(f>  +  S(f,)-f'(<l))]i6X)rdr 

JTc 

=  f  f  [/"((I  -  i)4>  +  K<f>  +  H))  -  f"W] dt (6<f))XT dr+ 

+  /  t  f”  ((1  -  +  *{<(>  +  dt  (s<f>){6X)T  dr 

mtf"'(s(l  —  t)(f)  -I-  st(<f)  +  8<f))  -f  (1  —  s)(f>)  dsdt  \8(i)^  XrdT 
_ 

+  /  /  f”{{1--i)<l>^-K<l>  +  ^<i>))dt(S<f>)(^X)Tdr  'i{8<t>,8X)£X. 

Jt  cJ  0 

Note  that  f{<f,)  =  03(6^^-^  -  /’((f,)  =  C3(Cie^i^  -b  /"(<!>)  = 

-  Cle-^^^)  and  /"'(c^)  =  -b  By  Lemma  2.2  we 

have  that  for  all  real  ntunber  m  >  1, 


ll/'(<?^)IU"*(rc)  < 


m^CUWl 
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ll/"Wllf(rc)<c{l  +  |r|  + 


m}cm\ 


in}' 


and 


lir(«lli~(rc)  <  c  {i  +  |r|  +  e  iri} "  , 


where  C  is  a  generic  constant  independent  of  (f>.  Trace  theorems  for  Q  C  implies 
that  for  all  p  >  1,  g  >  1  and  r  >  1, 

INIIinrc)  <  C||<^||i 

Il«||t.(rc)  <  v«€ir‘(n) 

and 

linii'cr.)  < Clkllw  Vr€ir»/^(r). 

We  fix  some  m>l^p>l,  q>l  and  r  >  1  with  ■^  +  p  +  ^  +  r~^’  Then  for 
every  (So,  ^ A)  €  X  we  have  that 

/  f  ((1  -  t)<f>  +  t((i>  +  S<f>))  dt  (S(I>)(6X)t  dr  < 

Jtc  •'0 

<  sup  II/" ((1  -  t)4>  +  t(<f>  +  ^fli))||L”*(rc)  ll<^</*llL'>(rc)  ll<^^llL»(rc)  Iklli.’-(rc) 

0<t<l 

<c{l  + in +  in}*  Mil  ll«||,  ||r||./j,r  Vr€ff'/^(r). 

Similarly,  we  have  that  for  every  S<f>  €  H^(ri)  and  every  r  € 

Iff  t  f" (s(l  —  t)(f>  +  st((f>  +  S<l>)  +  (I  —  s)<l>)  dsdt\6<f>\‘^  Xt  dr  < 

JTc  Jo  Jo 

<c{l+|r|  +  e’"’’^“‘|r|}"  M|?||A||i||T||,/3,r. 

Returning  to  (6.29)  we  obtain  that  for  all  r  €  and  (S(f>,  6X)  €  X, 

((2(?  +  «^,A  +  <A)-Q(.^,A),r)-  L  f'(<i>){6(f>)  XrdT-  L  /'(^)(iA)r<ir< 

<  C  {l  +  in  +  in}  ”  {||<^||.  ||«||.  +  11^,^111  ||A||i}  ||Tl|,/2,r , 


SO  that  we  conclude  Q  as  a  mapping  from  X  to  H  is  Frechet  differentiable 

and  its  derivative  Q'{<l>,  X)  is  given  by 

Q'(^,A)(^,<S),r)=^  [/"('^)^Ar  +  /'(^)iirl<ir  VT€^‘/"(r). 


Hence,  talcing  into  account  the  remarks  in  the  beginning  of  the  proof,  we  have 
justified  that  G  is  Frechet  differentiable  and  its  Frechet  derivative  G'((f>,  X)  is  defined 
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as  follows.  For  each  (^,  A)  €  X,  G'(^,  A)(^,d>)  =  (C^fj,6)  for  (^,a>)  €  X  and 
(C,  t},9)  eY  if  and  only  if 


(C.-)r  =  ^/  uj7rdr+  f  (f{<l>)^-arl;)ndr  V  tt  G  if'/2(r) , 

£>o  Jr  A  J^c 

{fj,u:)  = - f  rpivdO,  V  a;  € 

^  0  Jsi 


(6.30) 

(6.31) 


and 


{e,T)r=  f  /"(^)^Ardr+  /  (f\cl>)-a)u>Tdr  Vr€ir'/"(r);  (6.32) 

JTc  •'T'c 

or,  equi\'alently, 

c 


on  Ta‘, 

f'(<f>)ip  —  atl>  on  Tc] 
(  0  on  To  , 


fj  — - Ip  in  ft 

^0 


and 


1 


+  {f'{(f>)  -  a)u  on  Fc  ; 

0  onFoUryi. 


(6.33) 


(6.34) 


(6.35) 


(These  defining  equations  can  be  formally  derived  by  differentiating  (6.21)-(6.23).) 
It  is  easy  to  verify  from  the  above  equations  that  for  each  (^,t5)  €  X,  we  have 
(C)^)^)  €  Z,  i.e.,  G'(<^,A)  maps  X  into  Z;  furthermore,  using  trace  theorems  and 
Lemma  2.2  we  obtain  that 

IlCllo.r  <  +  ||/'(<;^)IU‘‘(rc)  II^IU4(r)  +  a||t^||o,r 

OO 


{i  +  in 


+  e 


c||<^|| 


'iri}’'"  W.+CMU, 


0  <  iillflli 


and 


IlCllo.r  <  Wf'WWiHTc)  II^IU«(r)  PllL«(r)  +  ll/'(<^)llL<(rc)  ll‘^llL‘‘(r)  +a||d;||o,r 


<  C  {l  +  |r|  +  e'""*"’  in}*'*  {Will  IIAII,  +  Hill, }  +  CIPIIi . 

Thus  G'[(f>,  A)  €  C{X-,  Z),  i.e.,  we  have  shown  that  (6.6)  hold. 

To  show  the  second  order  differentiability  of  G,  again  for  clarity  we  will  ex¬ 
amine  only  one  nonlinear  term  appearing  in  the  definition  of  G',  e.g.,  the  term 
r  i->  /"(<^)^At  dF.  We  define  a  mapping  R  :  X  C(X-,H~^^^(T))  by 
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f"{(f>)i,XT  dT  for  all  (<^,A),  (^,tD)  €  X  and  r  € 

For  each  given  (<^,  A)  €  X  we  have  that  for  all  {84>,  (5A)  6  X, 

( [i2((^  +  A  +  a)  -  i2(^,  A)](0,  a;),  r ) 

-/  /'"(^XWArdT-  /  /"(<^)^(a)rdr 

•/Fc  </rc 

JTc 

+  !  [n<i>  ^  s<f>)  -  r{<f>m8\)T  dr 

JTc 

mt  —  t)<f>  +  st(<f>  +  +  (1  —  s)(l>)ds  dtrj^\8(l>\^  XrdT 

_ 

+  /  /'  /"'((I  -  +  t{<i>  +  d<  {8(f>){8X)i>T  dr . 

Jt  c  ^ 

Thus  similar  to  the  analysis  ensuing  (6.29),  we  can  show  that  the  operator  R  is 
Frechet  differentiable  and  its  derivative  is  defined  by: 

i2UA)-((^,^),(t^,^))=  /  r(<f>)U^rdr+ f  r{<j>)xflTdr. 

JTc  JTc 

Hence,  G  is  second  order  Frechet  differentiable  and  G''{<i>,X)  is  defined  as  follows. 
For  each  (^,  A)  €  JT,  G"(4>,  A)  •  ((^ ,  w),  (0,5))  =  (<,  n,  0)  for  ((^,  a,),  5))  €  AT  x  AT 

and  (C,  fj.9)  eY  if  and  only  if 

(C,7r)r=/  f"(<f>)$rl^^dr 
JTc 

{fj, u;)=0  yujEH^iQ) 

and 

{9,r)r=  f  [f"\<f>)$^X  +  f"(<i>)^l]Tdr+[  f''{<f>)~^u>Tdr  VrG^'/2(r); 

JVn  JVn 


'Fc 

or,  equi%'alently. 


Fc 

'0  on  Tx ; 

C  =  '  ^  on  Tc ; 

,0  on  To , 

Tjf  =  0  in  fi 


Q^[  A  +  X]  +  Cb  on  Tc  ; 

0  on  To  U 


(6.36) 

(6.37) 

(6.38) 


and 


Furthermore,  using  Lemma  2.2,  (6.36)-(6.38)  and  trace  theorems,  we  may  derive  a 
bound  for  G"(4>,X)  for  each  given  (<f>,X): 


l|G"(^.  A)||y  <  C  {l  +  |r|  +  |r|}  (1  +  ||A||,) 

for  some  constant  C,  so  that  G"  is  bounded  on  every  bounded  subset  of  X.  I 

A  solution  (</>,  A)  of  the  problem  (5.2)-(5.3),  or  equivalently,  of  (6.27),  is  nonsin¬ 
gular  if  the  linear  system 


I  £7  grad  •  grad  dfl -f  /  f'(<f>)^rp 

Jsi  JTc 


dT 


+ 


Oo  Jta 


(6.39) 


and 


f  (jgradA  ■  gradwdfl -(-  /  f"(<l))  ^  Xu  dT  +  j  f'{(l>)Xujdr 
Jq  Jtc  c 

-  —  f  ^udQ  =  {Tj,u) 

^0  JQ 


(6.40) 


has  a  unique  solution  (<^,  A)  €  A"  for  every  Ci  ^7  €  if^(fi)*. 

An  antdogous  definition  holds  for  nonsingulax  solutions  of  the  discrete  optimality 
system  (6.2)-(6.3),  or  equivalently,  (6.28). 

It  is  e\*ident  that  (6.39)-(6.40)  has  a  tmique  solution  for  large  enough  <7,  e.g.^ 


a  >  max 


It  is  reasonable  to  assume  that  (6.39)-(6.40)  htis  a  unique  solution  generically  with 
respect  to  <r,  z.e.,  the  optimal  solutions  are  almost  always  nonsingular.  Thus  The¬ 
orem  6.1  and  Proposition  6.3  lead  to  the  following: 

Theorem  6.4  Assume  (^,  A)  is  a  nonsingular  solution  of  the  optimality  system 
(5.2)-(5.3).  Assume  that  the  finite  element  spaces  satisfy  the  condition  (6.1). 
Then,  there  exists  a  ^  >  0  and  /iq  >  0  such  that  for  h  <  ho,  there  exists  a  unique 
nonsigular  solution  {(f>^,X^)  of  the  discrete  optimality  system  (6.2)-(6.3)  satisfying 
—  <p\\i  -|-  ||A^  —  A||i  <  8.  Moreover, 

+  -^0  ash^O.  (6.41) 

If,  in  addition,  the  solution  of  the  optimality  system  satisfies  {<f,X)  €  if”*'*‘^(fl)  x 
ir"*+i(n)j  then  there  exists  a  constant  C,  independent  of  h,  such  that 


11^  -  +  ||A  -  A‘||,  <  Cfc^dl^ilU+i  +  ||A|U+,) .  I 


(6.42) 


A  consequence  of  Theorems  6.4  is  the  following  corollary  that  gives  error  esti¬ 
mates  for  the  approximation  of  the  controls. 

Corollary  6.5  Assume  (<f>,X)  is  a  nonsingular  solution  of  the  optimality  system 
(5.2)-(5.3).  Assume  that  the  finite  element  spaces  F*  satisfy  the  condition  (6.1). 
Define  the  approximate  control  by 

=  onTA. 

oo 

Then 

ll«‘ -  u||,/2,rx  ^  0  (6.43) 

If,  in  addition,  the  solution  of  the  optimality  system  satisfies  (^,  A)  €  x 

then  there  exists  a  constant  C,  independent  of  h,  such  that  for  h  <  ho, 

l|u‘  -  u||./2,r..  <  +  IIAIUh-O  .  (6.44) 

oo 

Proof:  Recall  that  u  =  “^A  on  Fa;  see  (5.1).  Then  (6.43)  and  (6.44)  follow 
trivially  from  (6.41)-(6.42)  and  the  inequalities  (see  [1]) 

11“  -  “^lli/2,rA  =  “  A'‘lli/2,rc  ^  ^11-^  “  '^*lli/2,r  <  ^l|A  -  A*||i ,  I 

Now  we  wish  to  apply  Theorem  6.2  to  derive  T^(rc)-error  estimates  for  the 
approximations  of  u.  To  this  end,  we  assume  the  domain  Q  is  convex  and  for  each 
given  e  6  (0, 1/4),  we  introduce  spaces 

H  =  X  and  W  =  x  . 

Note  that  X  C  H  with  a  compact  imbedding  so  that  (6.18)  implies 

||(T  -  T*)|U(y;H,  ^  0  asA^O. 

Again  using  finite  element  approximation  results  in  [6]  we  have  that  if  Q,  is  convex 
and  T{(:,r],d)  €  x  then 

l|(T-r‘)(C,7;.«)||H  <  CA"— 

Proposition  6.6  For  each  (^,  A)  €  W,  the  operator  G'{<j>,X)  :  X  Y  defined  by 
(6.33)-(6.35)  can  be  extended  as  a  linear  operator  of  C{H]Y).  Furthermore,  the 
mapping  w  — »•  G'{w)  is  continuous  from  W  onto  C{H',Y). 

Proof.  Note  that  W  C  L°°(Q,)  x  L°°{Q)  and  C  L^(F)  with  continuous 

imbeddings.  For  each  ((j>,  A)  €  W,  we  can  easily  verify  from  (6.33)-(6.35)  that 

||G'(<A,A)(^,a))||y 

<  ^(^){||w|io,rA  +  ll^llo,rc  +  ll’^llo  +  Iloilo, Fc} 

<  Ct  A(4,)  {||<5||i/2+,  +  |iiA||i/2+.}  V  (^,<i)  €  if , 


where 


A{(f))  =  C  max 


[-  - 
\  ^0  ’  Co  ’  |3c|< 


max  (l/'(x)|  +  a),  max 

ll«^ll3/2+«  |x|<|l^||3/2+ 


P||s/2+,)} 


and  Ct  is  a  constant  such  that  ||V’||o,r  ^  C'rj|t/’||i/2+€  all  ^  The 

desired  results  follow  easily  from  this  estimate.  I 

If  ((^,  A)  is  a  nonsingular  solution  of  (5.2)-(5.3),  using  the  denseness  of  H^(Q)  in 
and  regularity  theories  for  (6.39)-(6.40),  we  infer  that  (6.12)  holds. 

Thus  we  have  verified  all  the  requirements  in  Theorem  6.2  so  that  we  can  draw 
the  following  conclusion: 

Theorem  6.7  Assume  fi  is  convex  and  (<f>,  A)  is  a  nonsingular  solution  of  the 
optimality  system  (5.2)-(5.3).  Assume  that  the  finite  element  spaces  satisfy  the 
condition  (6.1).  Then,  there  exists  a  S  >  0  and  ho  >  0  such  that  for  h  <  ho, 
there  exists  a  unique  nonsigular  solution  (^^,  A^)  of  the  discrete  optimality  system 
(6.2)-(6.3)  satisfying  ||<^^  -  +  IjA^  -  Aj|i  <  6.  If,  in  addition,  the  solution  of 

the  optimality  system  satisfies  (<?i,A)  6  x  then  there  exists  a 

constant  C ,  independent  of  h,  such  that 

U  -  +  ||A  -  A'‘||,+a/2  <  C;i”*-‘+'/'“(||<A||m+i  +  IIAIIm+i) .  I  (6.45) 


A  consequence  of  Theorem  6.7  is  the  following  corollary  that  gives  the  L'^{Va)- 
error  estimates  for  the  the  approximation  of  the  controls. 

Corollary  6.8  Assume  Q.  is  convex  and  ((f>,  A)  is  a  nonsingular  solution  of  the 
optimality  system  (5.2)-(5.3).  Assume  that  the  finite  element  spaces  satisfy  the 
condition  (6.1).  Define  the  approximate  control  by 

=  onT^. 

Oo 

If  the  solution  of  the  optimality  system  satisfies  {(f>,X)  G  x 

then  for  each  e  €  (0, 1/4)  there  exists  a  constant  C ,  independent  of  h,  such  that  for 

h  <  ho, 

l|u‘  -  .xlloj.  <  +  l|A|U+i) .  (6.46) 

Oo 

Proof:  Recall  that  u  =  —j^X  on  Fa;  see  (5.1).  Then  (6.46)  follows  trivially  from 
(6.45)  and  the  trace  theorems  (see  [1]) 

P  -  “‘lloj..  =  r  11-^  -  ^  rll-'  -  ^  f  11-^  -  ■  ■ 

OO  ^0  Oo 
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